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Contributions of this Study 

Among the new significant aspects of the present work are (i) the 
treatment of the far-field boundary, (ii) the use of C-grid topology, with 
the branch-cut singularity treated analytically, (iii) evaluation of the 
effect of the envelope of prevailing initial states and, finally,- (iv) the 
ability to employ streakline/pathline 'visualization' to probe the unsteady 
features prevailing in vortex-dominated flows. The far-field boundary is 
placed at infinity, using appropriate grid stretching. This contributes to 
the accuracy of the solutions, but raised a number of important issues which 
needed to be resolved; this includes determining the equivalent 
time -dependent circulation for the pitching airfoil. A secondary counter- 
clockwise vortex erupts from within the boundary layer and immediately 
pinches off the energetic leading-edge shear layer which then, through 
hydrodynamic instability, rolls up into the dynamic stall vortex. The 
streakline/pathline visualization serves to provide information for insight 
into the physics of the unsteady separated flow. 


t This research is supported, in part, by AFOSR Grants (Nos. 87-0074 and 
90-0249), with supercomputer resources being provided by the Ohio 
Supercomputer Center. 
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Physical Characteristics and Background 


Specifically, the large amptitude rapid pitching motion associated with 
the initiation of a high angle - of - attack maneuver typically leads to the 
generation of a dynamic stall vortex whose evolution results in large 
transient lift, drag and moment that can, for short periods of time, produce 
loadings significantly larger than those expected during either steady, or 
quasi-steady flight. Indeed, the successful completion of an abrupt, 
drastic maneuver can depend upon the ability of holding the dynamic stall 
vortices in place, at least for the duration of the maneuver, and 
subsequently bleeding the excess accumulated vorticity in a controlled 
manner into the wake. Abrupt shedding of large amounts of locally 
concentrated vorticity can so rapidly alter the lift distribution on a body 
that a tumbling loss-of -control incident can occur, as the associated rapid 
changes In moment distribution cannot be tolerated. 

Recently, Carr (1988) has comprehensively reviewed the literature on 
the dynamic stall phenomenon and has also articulated the effect of key 
parameters on this phenomenon. Helin (1989) has also highlighted recent 
advances in the field, while stressing the importance of unsteady 
aerodynamics for highly maneuverable and agile aircraft. In addition, he 
has raised the important issue of the effect of flow separation on the 
formation of the energetic dynamic- stall vortex. These two reviews 
adequately point out some of the unresolved issues associated with the 
problem of dynamic stall. 

On the Analysis of Dynamic Stall 

The unsteady Navier- Stokes analysis of K. Ghia, Osswald and Ghia (1985) 
and Osswald, K. Ghia and U. Ghia (1986) Is modified to permit arbitrary 
three degree-of -freedom maneuvers, using body-fixed coordinates and a C-grid 
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topology. This formulation not only permits pitching motions, but also 
plunging and in-plane accelerating or decelerating motions; typically, the 
airfoil is pitched about the quarter- chord axis. The problem is formulated 
using vorticity and stream function as dependent variables in a non-using 
body-fixed reference frame in generalized coordinates. This formulation 
offers the important advantage over the primitive -variable formulation that 
the form of the governing equations in inertial and non- inertial reference 
frames is identical. The far-field boundaries are located at true infinity 
for this subsonic flow with its fully elliptic nature. The conformal 
mapping techniques used lead to analytical determination of the 
corresponding inviscid flow; this inviscid flow constitutes the true far- 
field boundary condition; by contrast, the studies of Visbal and Shang 
(1988) and Ekaterinaris (1989) place the far-field boundary at a finite 
distance from the airfoil, and employ free-stream conditions on the upstream 
boundary and zero streamwise gradients downstream. The present study uses 
an analytically determined clustered conformal grid, thereby avoiding 
numerical error in the computation of the metrics. The C-grid topology 
employed introduces a singularity at the trailing edge (TE) and all along 
the branch cut. For the latter, this singularity is treated using the 
method of analytic continuation, as developed by Osswald, K. Ghia and U. 
Ghia (1985) . The conditions of zero slip and zero normal velocity at the 
surface of the airfoil are implemented appropriately in terms of the stream 
function and vorticity. At the TE, the singularity in the grid does play a 
role in determination of the stream function, which is obtained by 
satisfying the Kutta condition. The vorticity at the TE is determined using 
the analysis of Osswald, K. Ghia and U, Ghia (1989). The direct numerical 
simulation (DNS) methodology developed by the authors is used to solve the 
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vorticity transport and stream function equations. Central differences are 
used for all spatial derivatives and no artificial dissipation is added 
explicitly . 

Results and Discussion 

In the present study, simulations are carried out for a NACA 0015 

airfoil undergoing constant 0-pitch-up motion. Two flow configurations are 

3 A 

attempted and they are for Re = 10 and 10 . Configuration I, with the 

3 

lower Re ~ 10 , is used in the development phase, since it permits the use 
of a smaller grid (274, 76) of which 114 grid points are placed on the body. 

On the other hand, configuration II, with Re = 10^ , is used to compare the 
results of the experiments of Walker, Helin and Strickland (1985) who 
considered Re * 45000. This latter configuration was run using a (444,101) 
grid with 204 grid points on the body; the size of the grid was selected 
based on the results of Visbal and Shang (1988) who had carried out a grid 

study and selected this size. In addition, the same constant 9-pitch~up 
motion as used by them is also implemented here and corresponds to 
nondimensional pitch rate 0.2 with nondimens ional acceleration time 

t Q “ 0.5, and pivot axis location measured from airfoil leading edge 
x 0 - 0.25. 

Results of configuration I, in Fig. 1, show that the dynamic stall 
vortex with its clockwise spinning fluid evolves as the shear layer from the 
leading edge is pushed away by the counterclockwise spinning vortex close to 
the body surface and subsequently the shear layer rolls up and forms a 
dynamic stall vortex. It appears that the eruption of counterclockwise 
spinning vortex from within the boundary layer is important to formation of 
the dynamic stall vortex near the leading edge. 
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Results for configuration II, in Fig. 2, show that 9-0 run exhibits 
unsteady results as opposed to the steady-state results obtained by Visbal 
and Shang. Their steady-state solutions could strictly be a consequence of 
the use of both explicit as well as implicit smoothing, i.e., artificial 
viscosity, to maintain stability in their numerical calculations, which 
employed the method of Beam and Warming. Present results for configuration 
II, although not depicted here, show that there are grid related 
oscillations near the leading edge in the vorticity contours and grid 
structure and perhaps its density needs to be altered before generating new 
results and analyzing them. Still, in this case, also the secondary 
counterclockwise vortex erupts from within the boundary layer on the surface 
to form the dynamic stall vortex. Unlike wind-tunnel tests, the numerical, 
experiment in the present study computes the vorticity field directly. By 
evaluating various individual terms in the vorticity- transport equation, it 
is possible to examine vorticity accumulation and generation at the body 
surface as well as in the flux from the boundary to reveal the underlying 
mechanism and the role of unsteady separation on the evolution of the stall 
vortex. This is possible once a comprehensive set of results are obtained. 

In summarizing, the constant 0-pitch-up experiment of Walker et al. 
(1985) is simulated using direct numerical simulation and an unsteady NS 
analysis. The preliminary results obtained so far provide the flow 
structure and the evidence that eruption of secondary counterclockwise 
vortex near the quarter chord point triggers the formation of the dynamic 
stall vortex. However, additional results are essential to obtain the 
budget of vorticity dynamics and to shed further insight into this mechanism 
underlying the evolution of dynamic stall vortex just stated and its 
relation to unsteady separation. Based on the existing results of a video 
presentation of the numerically simulated evolution, convection and shedding 
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of the dynamic stall vortex is created, but for quantitative information 
comprehensive data is needed. It should be added that, in the earlier 
results for flow past a static Joukowski airfoil at 53° , the authors have 
seen the formation of a secondary counterclockwise vortex before the leading 
edge shear layer forms a large clockwise spinning vortex. 
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RING 



FIG. 1 CONT. CONTOURS OF INSTANTANEOUS STREAM FUNCTION AND VORTICITY DEPICTING DYNAMIC STALL 
Re = 1,000, (a-b) a = 40.007°, (c-d) a = 32.295. 
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FIG 2 VORTICITY CONTOURS FOR FLOW PAST A NACA0015 AIRFOIL, 
Re = 10,000, a = 0°, t = 8.00. 
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FIG. 2 CONT. ENLARGED VIEW OF LEADING EDGE DEPICTING EVOLUTION OF 
DYNAMIC STALL, Re = 10,000, a = 13.651°. 
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INERTIAL VS BODY FIXED OBSERVER 



Inertial Observer 


Position of Fluid Particle: 
Velocity of Fluid Particle: 
Acceleration: 

Vorticity: 


n 

V, 


a i 

ui = V x Vi 


Body Fixed Observer 
(Apparent) 

f 

V 

a 

« = vxy 


KINEMATICS 

?i = r B /i(t) + f 
Vi = V B /j{t) + V + X r 

aj = a B /i(t) + a + 2 ti B (t) x V + a B (t) x f + fts(t) x ^b(*) x f 

u>r = V x { V B /i(t ) + V + x f} 

= 0 + VxV + 2fi fl (l) 

— uj -J* 


Arbitrary Maneuver Defined by: r B /i{i) 

v*//W = ^ MO 

« b //(0 = MO = 
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UNSTEADY INCOMPRESSIBLE 
NAVIER- STOKES EQUATIONS 

INERTIAL OBSERVER 
PRIMITIVE VARIABLES 

Continuity 

V • Vf = 0 

Linear Momentum (Bernoulli’s Form) 

^+(VxV,) X V,+ ^(VxVxVi)=-v(p + L|^) 
VELOCITY- VORTICITY 

Continuity 

V • Vi = 0 

Kinematic Definition of Vorticity 

V X Vj = L)[ 

Vorticity Transport 

^ + Vx(w / xi’,)+^(VxVx« / ) = 0 
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UNSTEADY INCOMPRESSIBLE 
NAVIER-STOKES EQUATIONS 


BODY FIXED (APPARENT) OBSERVER 
PRIMITIVE VARIABLES 

Continuity 

V • V = 0 

Linear Momentum (Bernoulli’s Form) 

^ + a B (t) xf+(Vx P) x V + 2 (V x V x v) = 


—V p + 


V 


+ a B /i(t ) • r - 


B(t ) x r 


VELOCITY- VORTICITY 

Continuity 

V • Vi = 0 

Kinematic Definition of Vorticity 

VxV/=u/ 

Vorticity Transport 

^ + V x (u>/ x V) + ^(V x V x w,) = 0 

Kinematic Relationship Between Apparent and Inertial Velocity 
V = V, - V B/ i{t) - ft B (t) x f 

• Solve For Inertial Velocity and Inertial Vorticity 
Directly in Body Fixed Frame 

• Form of Governing Eqs. Unaltered 


140 


Inertial Velocity Boundary Conditions Remain Unaltered 


Only Differences Are 

Inertial Vorticity Advects with Apparent Velocity 

Additional Vorticity is Created at Body Surface Due 
to Acceleration of Body 



For Two-Dimensional Flow, use 


DISTURBANCE STREAM FUNCTION- VORTICITY 
FORMULATION IN GENERALIZED BODY FIXED 
COORDINATES 


Definition of Disturbance Stream Function (Deviation from Uniform 
Flow) 


0 INERTIAL = [y + ®o] + 


Where 'Fo is at yet an Unknown Integration Constant Representing 
a Displacement of the ZERO STREAMLINE at INFINITY 


A 

<1 







Arbitrarily Maneuvering 

iPinertial = [x 2 cos $(t) - a : 1 sin 0(t) + ’Fo] + *) 
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ELLIPTIC STREAM FUNCTION PROBLEM 


1 x 1 /«'• «f' 5 \ 

3f‘ \y/s W ) d? \jg~i ~ ) ~ ~' / ° U ' 


Subject to the Boundary Conditions 

ip? rs = 0 at INFINITY 

4>f' s = {i 2 [V'i / ,(()-cos«(I)]-i , [V'i / ,(()-sin«(;)]-^((xi) 2 + (i I ) J ]} 

4 

Along Body Surface 


V>/. = y 

P 

0 - 

* 

Vo 


75 = o 

at INFINITY 



STRICTLY AN INVISCID EFFECT (ALL VISCOUS 
DISTURBANCES DISSIPATE WELL BEFORE INFINITY) 

DIRECTLY REPRESENTS UNDERLYING INVISCID 
CIRCULATION SET BY AN INVISCID KUTTA CONDITION AT 
TRAILING EDGE 
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KUTTA CONDITION FOR <!>„ 
(WEDGE TRAILING EDGE) 



lim 

~TE 

BRANCH CUT 


{ 


1 d^jj^yjscip 

s/gl 7 d? 


} = 


[{cos 0(t) - Vgp(t) + x 2 n B (t)}e 1 + (sin0(<) - V^j(t) 


Branch 

Cut 


— x 1 n B (<)}e 2 • t 
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VORTICITY TRANSPORT 


^tr + ^ H(vsn) + ^ {^Mv 2 }) = 

— S— ( , _d_ ( ffn d<^/\ 


where 

^,1 [{cos g(<) - V£ fI (t) + x 3 n B (t)}ii + {sin OJt) - %,(<) - adfl B ( <)}e 2 ] ■ c t ^f /s 

Wv^J + <9£ 2 

y /j r a _ [{ cosg (0 ~ ^//(O + » 2 «/»(0}e'i + {sinfl(Q - !*„(<) - ^O^)}^] • e 2 q^DIS 

922 /V 9 + ac 


Subject to the Boundary Condition 

<jj] = 0 at Infinity 


and Along the Body 


d 


—\/gu)r = — 
vy 1 


(gndW S _) . 3 (g u ty?' s 

\Vs 80 ) ae Vv5 30 


) 


Subject to the Constraint 


ffn dtf'S 

V9 9? 


[{cos^(O-V , i // (t)+x 2 O B (O}e 1 +{sin0(<)-yj // (t)-x 1 fis(O}e2]-e 2 
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THE PITCH UP MANEUVER OF 
VISBAL AND SHANG 


A 

d / 



= 0 
= 0 

L 


o B (t) = n 0 [*H^)(i-exp-(^)o] 

MO = -fy, (l — exp — (-y— 

... . 4.6 . 4 . 6 . 

MO = -(— )fi 0 exp-(— )f 


to 


Where 

n 0 - Nondimensional Pitch Rate; 0.2 

t a - Nondimensional Acceleration Time; 0.5 

x 0 - Pivot Axis Location Measured from Airfoil LE; 0.25 
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